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o ■ 
o : 

CN ■ Abstract 
S_i ' 

Q H | We derive the late-time behaviour of tilted Bianchi VIIo cosmologies with an irro- 

■ tational radiation fluid as source, and give the asymptotic form of the general solu- 
^-j- . tion as t — > +oo, making comparisons with the dust-filled models. At first sight the 

I | radiation-filled models appear to approximate the flat FL model at late times, since the 

■ Hubble-normalized shear and the tilt tend to zero and the density parameter tends to 
one. The Hubble-normalized Weyl curvature diverges, however, indicating that phys- 
ically significant anisotropy remains. We also discuss the influence of a cosmological 

■ constant on this phenomenon. 

o : 

1 Introduction 

The evolution of spatially homogeneous cosmologies of Bianchi type VIIo with a non-tilted 
perfect fluid source and zero cosmological constant has recently been shown to exhibit a 
number of new and interesting features at late times (see Wainwright et al 1999 and Nilsson 
et al 2000) 1 Firstly, they exhibit the breaking of asymptotic self- similarity, and thus serve 
as counter-examples to the hope that any Bianchi universe is approximated at late times by 
^ | an exact self-similar Bianchi universe (see WHU, pages 2578-2579). 

The breaking of asymptotic self-similarity is characterized by oscillations in the Hubble- 
normalized shear that become increasingly rapid in terms of the cosmological clock time t as 
t — > +oo. Depending on the equation of state, this behaviour also leads to Weyl curvature 
dominance, which refers to the fact that Hubble-normalized scalars formed from the Weyl 
curvature tensor become unbounded as t — > +oo (see WHU, page 2588). 

Asymptotic self-similarity breaking also has interesting implications as regards isotropiza- 
tion. In a spatially homogeneous cosmology, there are two physical manifestations of anisotropy: 



(a) the shear of the timelike congruence that represents the large-scale distribution of the 
matter; 

(b) the Weyl curvature, which can be viewed as describing the intrinsic anisotropy in the 
gravitational field: it determines up to four preferred directions, the principal null 
directions of the gravitational field. 



1 From now on we will refer to these papers as WHU and NHW respectively. 
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The dynamical significance of these anisotropies can be quantified by normalizing the shear 
tensor and the Weyl curvature tensor with the Hubble scalar, and defining two scalars, the 
shear parameter E and the Weyl parameters Wi, W2 (see Section 0] for the definitions). 2 

Bianchi I cosmologies isotropize for all values of the equation of state parameter 7 — 1 = 
p/p, i.e. < 7 < 2, in the sense that E — > and Wi,W2 — > as t — ► +00. Indeed 
these models are asymptotically self-similar, and the asymptotic state is the flat Friedmann- 
Lemaitre (FL) model. On the other hand, for the class of Bianchi VIIo cosmologies the value 
7 = |, corresponding to a radiation fluid, acts as a bifurcation value as regards isotropization: 
if 7 < |, the shear parameter E^Oast— > +00, while if ~ < 7 < 2, E oscillates indefinitely 
away from zero, and does not approach a limit (see WHU, page 2581). In addition, the 
isotropization that occurs in the case 1 < 7 < | differs significantly from the isotropization 
that occurs in Bianchi I cosmologies, due to the previously mentioned phenomenon of Weyl 
curvature dominance. Specifically, in Bianchi VIIo cosmologies, if 7 = 1, the Weyl scalars 
Wi and W2 are bounded but do not tend to zero as t — ► +00, while if 1 < 7 < |, Wi and 
W2 are unbounded as t — > +00 (see WHU, page 2581). We say that shear isotropization 
occurs if 7 < 3, but that Weyl isotropization does not. 

More recently, Coley & Hervik 2005 have investigated the effect of tilt on the dynamics 
of Bianchi VIIo cosmologies at late times. In general tilted Bianchi cosmologies have three 
additional degrees of freedom which describe the orientation of the fluid four-velocity relative 
to the normal to the group orbits (see King & Ellis 1973). For simplicity Coley & Hervik 
restricted their attention to models with one tilt degree of freedom. They found that the 
value 7 = I of the equation of state parameter acts as a bifurcation value as regards the 
tilt: if 7 < |, the tilt variable tends to zero as t — * +00 and the model undergoes shear 
isotropization as in the non-tilted case, while if | < 7 < 2, the tilt variable does not tend to 
zero as t — > +00. The analysis of Coley & Hervik does not include models with 7 = |, the 
bifurcation value. Our first goal in this paper is to give a proof of the asymptotic behaviour 
of this class of cosmologies, i.e. the radiation-filled Bianchi VIIo cosmologies with one tilt 
degree of freedom. Our second goal is to give the asymptotic decay rates as t — > +00 for the 
various physical quantities. 

The outline of the paper is as follows. In Sectional we present the evolution equations for 
Bianchi VIIo cosmologies with a tilted irrotational perfect fluid. In Section |3] we establish 
the limits of the dimensionless gravitational field and matter variables at late times for the 
case of a radiation fluid with zero cosmological constant. Section 0] concerns the question of 
isotropization at late times. Section|5]presents the asymptotic form of the solutions. Sectional 
considers the effect of a positive cosmological constant on the results of the previous sections. 
We conclude in Section [7| with a discussion of the implications of the results. The details of 
the proof for the limits are given in Appendix El The derivation of the asymptotic forms is 
given in Appendix [BJ 

2 The evolution equations 

As in previous recent investigations of Bianchi VIIo cosmologies (see WHU, NHW, Coley 
and Hervik 2005), we write the Einstein field equations using Hubble- normalized variables 
within the framework of the orthonormal frame formalism of Ellis and MacCallum 1969, in 

2 We shall use Wi and W2 in this paper instead of W as in WHU and NWH because Wi and W2 are 
spacetime scalars while W is defined with respect to the congruence, and hence is not a spacetime scalar. 
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which the commutation functions Yab of the orthonormal frame {e a }, defined by 



[^a, ®i>] T ab^c , 

act as the basic variables of the gravitational field. 3 We choose the frame to be invariant 
under the isometry group G3, and eo = n, the unit normal to the group orbits. It follows 
that 7 c a b depends only on t, the clock time along eo- The non-zero components of the 
commutation functions are (see Wainwright and Ellis, page 39) 

{H, Cap, n aP, da, 

The variable H is the Hubble scalar, a a p is the rate of shear tensor, n a p and a a describe the 
spatial curvature, and fl a is the local angular velocity of the spatial frame {e a } with respect 
to a Fermi-propagated spatial frame. The energy-momentum tensor in the case of a tilted 
perfect fluid is 

T ab = pu a u b + p(g ab + u a u b ) , 

and we assume linear equation of state p — (7 — l)p. The fluid 4- velocity u is written in the 
form 

V 1 — VbV° 

where the spacelike vector v is orthogonal to n, and satisfies < v b v b < 1. The vector v 
is called the tilt vector of the fluid. It has components (0, v a ) relative to the orthonormal 
frame {n, e a }. Writing the energy-momentum tensor in the form 

Tab = pn a n b + 2q (a n b ) + p(g a b + n a n b ) + vr ab , 

the source terms (p,p,q a ,^ap) can be expressed in terms of p and v a (see Appendix A of 
Hewitt et al 2001). The Hubble-normalized commutation functions and the density paramter 
VL are defined by 

(E Q/3 ,iV a/ j, A a ,R a ) = {o- af3 ,n al 3,a a ,rL a )/H , tt = p/(3H 2 ). (1) 
We introduce the usual dimensionless time variable r according to 

- = -■ (2) 

For more details, see Appendix A of Hewitt et al 2001. 

For the class of Bianchi VII cosmologies with a tilted perfect fluid, we have A a = 0, and 
one of the eigenvalues of N a p is zero. We choose the spatial frame {e Q } so that ei is the 
corresponding eigenvector. It then follows that 

N la = , (3) 

and the evolution equations of the Ni a (Hewitt et al 2001, equation (A. 12)) imply 

R 2 = E 13 , R 3 = -S12 • (4) 

3 Indices are with respect to the basis of the frame vectors instead of the coordinate basis; lower case Latin 
indices run from to 3; Greek indices from 1 to 3; upper case Latin indices from 2 to 3. 
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We next impose the restriction that the perfect fluid is irrotational, which is equivalent to 4 
setting 

V2 = v 3 = = E 12 = S13 • (5) 
We now label the shear variables E a/ 3 and the spatial curvature variables N a p as follows 

\ /0 

;S^)=| S + + v/3S- \/3£ x , (N aP )= iV + + v^AL v^AL 




v^Sx S+-v / 3S_/ \0 v^AL N+-V3N-, 

(6) 

and also drop the index on v\ and 

With these restrictions, the general evolution equations (A.11)-(A.13) and (A. 29) in 
Hewitt et al 2001 reduce to the following system of differential equations: 

E' + = (g-2)S + -2(A^! + ^)-G; 1 7^ 2 (7) 

E'_ = (g-2)E_-2i?E x -2A+AL (8) 

E' x = (g-2)E x + 2i?E_ - 2A+A X (9) 

A^ = (g + 2E + )A + + 6E_A_ + 6E X A x (10) 

Ai = (g + 2E+)A_ - 2i?A x + 2A+E_ (11) 

A^ = (g + 2E + )A x + 2iiAL + 2A + E x (12) 

v' = Gl 1 [(3r - 4) + 2E+] (1 - «> , (13) 

where 

G± = 1 ± (7 - l)v 2 . 
The density parameter Q is given by 

n = 1 - £+ - £ 2 - E 2 X - A! - A 2 > , (14) 

and the deceleration parameter q is given by 

g = 2(E^ + E 2 _ + E 2 ) + ^[(37 - 2) + (2 - 7 )<; 2 ]fi . (15) 

There is one constraint left from (A. 27) of Hewitt et al 2001: 

= 2(E_A X - E x AL) + G'^Qv . (16) 

The evolution equation (A. 28) of Hewitt et al 2001 for Q simplifies to 

Q' = [2q - G; 1 [(3 7 - 2) + (2 - 7 )w 2 ] + 2G- 1 ^ + v 2 } Q . (17) 

The variable R in equations (|7jl— (|13p. which describes the rate of which the spatial frame 
is rotating in the 23-space, is unrestricted. We will eliminate it by introducing rotation- 
invariant variables as in (j20j) - (f2*2*J) below. 

The inequality (fT4*j) implies that the variables E + , E_, E x , A_ and A x are bounded. 
The remaining gravitational field variable, A + , however, is not bounded. In fact, it has been 

4 This result follows from the expressions for the fluid vorticity W a in Hewitt et al 2001, Appendix C. See 
King & Ellis 1973 for the original result. 
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shown (see equation (4.2) of Coley and Hervik 2005) that 



lim N + = oo 

T— >+00 



for all initial states with Q > and | < 7 < 2. We thus introduce a new variable 

M=±. (19) 

It can be seen from equations (jHJ) & (fTTj) . and & (|T2"|) that the unboundedness of iV + leads 
to increasingly rapid oscillations in the pairs of variables (£_, N_) and (S x , N x ). We thus 
introduce an angular variable ip and two rotation-invariant variables TZi and 7?.2 as follows 

Tlx = E 2 . + £ 2 + iV 2 + iV 2 (20) 
ft 2 cos 2^ = Y?_ + £ 2 - Nl - Nl (21) 
^ 2 sin2^ = 2(S_iV_ + S X A^ X ) . (22) 

It follows from (J20D (E2D that 

ft 2 - ft 2 = 4(£_iV x -S x iV_) 2 . (23) 

The resulting evolution equations for the variables (S + , 7^i, 7^. 2 , M, if), v) do not contain the 
variable R. 

We note that the conditions 

S x =0, iV x =0, ,R=0, (24) 

which imply v = on account of ([16)1 . describe the non-tilted models. With these restrictions 
equations (|7)l- (|T3|l reduce to (3.7), (3.9) and (3.10) in WHU. In addition the angular variable 
ip in (|21j) and (|22j) coincides with the angular variable ip defined by equation (3.15) in WHU. 
This conclusion follows from noting that implies TZi = 1Z 2 and that 1Zi = -R^hu- 

3 Limits for the radiation-filled case at late times 

In this section we determine the limits of the variables 

(S + , Hx, TZ 2 , M, iP, v) . (25) 
for the case of a radiation fluid, given by 

7 = I • (26) 
For this value of 7, the expression (J15)) for q simplifies to 

q = 1 + X 2 + TZ 2 cos 2i\) , (27) 
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on eliminating Q using (|14|). In terms of the variables f)25j) the evolution equations (fTj)- (fT3*|) 

assume the form 5 



= -(l-£^)£ + -72i + (l + £ + )72 2 cos2^--^^ (2* 

n[ = 2[(1 + £+)£+ + 72 2 cos 2^] - 2(1 + £ + )(cos2^)72 2 (29) 

7e' 2 = 2[(1 + £+)£+ + TZ 2 cos 2^] 72 2 - 2(1 + £+)(cos 2^)72i (30) 

M' = -[(l + S + ) 2 + ^ 2 cos2^ + 3^ 2 Msin2^]M (31) 

^ = A + (l + £ + )fi s in2^ (32) 

v >=*Y, + (l-v 2 )v. (33) 



3-u 

Using (|2*3|) . the constraint (fT^jl becomes 



Equation (|14jl for the density parameter assumes the form 

Q = 1 - £+ - Tlx . (35) 
The evolution equation (fTTjl for Q becomes 



4?> 



2 



S;+^ 2 COs2V>+- — 

3 + 



. (36) 



Observe that equations (|2*5Jl - (j3l)|) are invariant under the interchange v — > —v, and that 
v = is an invariant set. We can thus without loss of generality assume that v > for tilted 
models. The above formulation of the evolution equations has one unsatisfactory feature, 
namely the presence of the term TZi/TZ 2 in the evolution equation (J32|) for if). By (|20|). 
72-1 > for all r, and on account of (|34|). 72i = for some r is not possible since then 
f2t> = 0, which is excluded for tilted matter-filled models. So TZ\ > for all r, i.e. 72-1 = 
is an invariant set, in the boundary of the physical state space. On the other hand, 72-2 can 
become zero, in which case the change of variables (|21j) - ([2"2"|) is singular and the evolution 
equations (J33J1 break down. Numerical simulations using the evolution equations (jJJ- 
()13|) . however, provide evidence that for typical initial conditions TZ 2 can become zero for at 
most a finite number of times. In other words for a typical choice of initial condition, there 
exists a time r such that 1Z 2 (t) > or TZ 2 (t) < for all r > r . Note that the evolution 
equations (|2*5|) -(|3^ |) are invariant under the transformation ip — > ip + | and 72-2 — > — 7Z 2 . 
We make use of this symmetry to assume, without loss of generality, that TZ 2 (r) > for all 
t > tq. In other words, for typical initial conditions we expect that for r sufficiently large, 
the transformation (|20|) - (|2*2*|) will be non-singular, and the evolution equations (j2*Hj) - (|3*3*|) will 
be well-defined. 



These equations can also be obtained from the evolution equations (4.10)-(4.18) in Coley and Hervik 
4 
3 ' 



2005, by setting 7 = 4, and noting the following relation between our variables and theirs 
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The limit of M is 

lim M = , (37) 

which follows immediately from (jl8j) . Since i/j' in (|32j) is dominated by M _1 , the trigonomet- 
ric functions in (|28j) - (J3H|) will oscillate increasingly rapidly as r — > +00. In order to control 
these oscillations, we modify the variables S + , T^i, M, Z and (in analogy with equations 
(17)-(21) in NHW) according to 

£+ = £+ - |M(1 + S + )^ 2 sin2^ (38) 

^1 = 1 ; — \~t? (39) 

1 + |M(^! - l-E + )^sin2^ 

M 

M = = (40) 

1 - ±M7e 2 sin2V> v ; 

= . (41) 

1 i 1 1\/TT} ™„ o„/, v 7 



2 



M^ 2 sin 2^ 



It follows from (|2*5|l . (|2*5jl. (}3~T]) . and (|3l)|) that the evolution equations for the new 

variables have the form 6 



e' = -n, - s + (i - £5.) - _ _ (1 - 25. - + M5 S+ (42) 

3 + 1> z 

K' x = 2 [£+ + + MB^J ^1 (43) 

M' — — [(1 + £+) 2 + M5 A? ] M (44) 



3 + v 1 



tt , (45) 



where the B's are bounded as r — > +00. Equations (J42|) - (j45j) depend explicitly on ijj through 
the B's which are all multiplied by a factor M. This fact, in conjunction with ()37j) . implies 
that the oscillatory terms in tend to zero as r — > +00. Observe that the variable 

7I2 only appears in the 5-terms, which means that it is not necessary to introduce an 
evolution equation for a modified variable TZ2- The evolution equation (}3*3*j) for v remains 
unchanged. 

The first step in the analysis is to use the evolution equations ()28|) and (|33[) for S + and 
v to prove that lim T _ >+00 u = 0. The structure of then allows the limits of S + and 

IZi as t — ► +00 to be determined, and subsequently the corresponding limits of £+, 1Z\ and 
72-2 • The main result is contained in the following theorem. 

Theorem 1 Any solution of the system i2fy) - \34\) with Q > 0, 1Z\ > and IZ2 > satisfies 

lim (S + , n u n 2 , v, M) = (0, 0, 0, 0, 0), (46) 

T— >+00 

with 

7Z 7Z 



6 We have chosen to write equation l|45|l with S + unbarred in the second term, since this form is required 
in the proof in Appendix 1X1 The factor of IZi in front of is crucial, and IZ2 has to be written as (lj^ L )lZx 



to produce this factor. We recall from (j34(l that — 1 < < 1. 
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Proof. The details of the proof are given in Appendix [XJ D 



4 Isotropization of radiation-filled models at late times 

In this section we will consider the physical implications of Theorem Q as regards isotropiza- 
tion of the radiation-filled cosmological models with zero cosmological constant. As men- 
tioned in the introduction, the anisotropy can be quantified by the shear parameter E and 

a/3 

the Weyl parameters. The shear parameter is defined by £ 2 = °^ H ^ , or equivalently, on 
account of ((H), by 

£ 2 = |£ a/3 E^ . (48) 

It follows from © and (@BJ) that 

E 2 = T? + + E 2 + E 2 . (49) 
We define the Weyl parameters Wi and W2 by 

(~i riabcd r~i * r^abcd 

Wi = W9 = bcd (50) 

1 48# 4 ' 2 48# 4 ' { } 

or equivalently, 

E a fjE — H a pH E a pH 
1 6# 4 ' 2 3# 4 ' 1 J 

where E a p and H a p are respectively the electric and magnetic parts of the Weyl curvature 
tensor with respect to the congruence orthogonal to the G3 group orbits (Wainwright and 
Ellis 1997, page 19). We define the Hubble-normalized components by 

r _ E a p _ H a p 

o a p - , n a p - -jp- , [oz) 

labelling them as with E a ^ in ©. It follows that 

Wi = S 2 + + £ 2 + £l - H 2 + - Hi - Hi (53) 
W 2 = 2{£ + H + + EJH- + EyUy) . (54) 

They are given by 7 

£+ = {! + £+)£+ - (E 2 _ + E 2 ) + 2(N 2 + iV 2 ) + ^G^tty 2 (55) 

£_ = {!- 2E + )E_ + 2N+N_ (56) 

£ x = (1-2E + )E X + 2N + N X (57) 

H+ = -3(£_AL + E x iV x ) (58) 

H- = -3E+AL - 2iV + E_ (59) 

H x = -3E + iV x - 2iV + E x . (60) 



7 See Wainwright and Ellis 1997, page 35 for the general expression for E a p and H a p. See Lim 2004, page 
49 for the expression for the individual components, but note that the components are normalized by 3H 2 
there. 
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The limits of the anisotropy scalars S, Wi and W2 as r — > +00 for tilted irrotational 
radiation-filled Bianchi VIIo cosmologies are given in the following theorem. 

Theorem 2 For any tilted irrotational radiation-filled Bianchi VIIq cosmology with zero 
cosmological constant, satisfying f2 > 0, > ; 7?.2 > and v > 0, the density parameter 
satisfies 

lim tt = l , (61) 



the tilt variable v satisfies 

lim v = , (62) 



T— > + 0O 



and t/ie anisotropy scalars satisfy 



lim ._S = 0, (Wi,W 2 ) = -^|[(cos2^,sin2^) + 0(M)] , (63) 



with 

lir ? TH = + °° • ( 64 ) 

Proof. It follows from flU, and (gHJ) that 

S 2 = + I (^1 + ft 2 cos 2V>) • (65) 

Equations (|61|1. (jB^j) and (p)3*j) now follow directly from Theorem ^ and equation (|33j) for 
Q. The asymptotic forms for Wi and W2 follow directly from (|53|) - (}60|) and the definitions 
(JIHJ)-(|221) of M, TZ\ and 7Z2, on noting that the terms containing N + are the dominant ones. 
The limit (164)1 is part of Theorem ^ □ 

The limits in Theorem El show that tilted irrotational radiation-filled Bianchi VIIq cos- 
mologies isotropize with respect to the shear, but not with respect to the Weyl curvature. 



5 Asymptotic solution at late times 

In this section we derive the asymptotic form as r — > +00 of the basic variables (S + , TZi, 7I2, M, v), 
which determine the asymptotic form of the anisotropy scalars. We begin by considering 
the so-called truncated differential equations for the variables S + , TZi and v, obtained from 
PS]) . (jjHj) and (J33|) by setting to zero the terms that involve M. We also write TZi = R 2 , 
since this facilitates the analysis. The resulting differential equations are 

K = -t + + ±l-R 2 - -±^{1 -±l- R 2 )v 2 (66) 

3+ir 

R' = £ + (l + ± + )R (67) 

v' = ^-^(1 - v 2 )t + v , (68) 
3 — v l 

where we use a hat to denote the solutions of the truncated differential equations. In Ap- 
pendix |B] we show that the solutions of (JHEJ)-((HBJ) exhibit a power law decay to zero. Since 
M decays to zero exponentially fast (see (jl49|0 it is plausible that the solutions £+, TZ\ and 
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v of the full system (J4*2~j) . (jlSj) and (jHHj) will have the same asymptotic form as the solutions 
of the truncated differential equations ([66 | -()68j ) . Finally (J3S|) and fj39j) show that S + and 
IZi have the same decay rates as E + and TZ\. An asymptotic expansion for M can most 
easily be obtained by using the auxiliary variable Z, defined by (j!22|) . with Q given by ()35|) . 
Knowing the asymptotic expansions for v and £+, and that M decays exponentially (see 
()149|) ). we can integrate ()125|) to obtain an asymptotic expansion for Z and hence recover Z 
from (|124|) . This process introduces another arbitrary constant Cm- 

These considerations lead to the following asymptotic form of the basic variables in any 
tilted Bianchi VII cosmology with irrotational radiation fluid. We give the asymptotic 
expansion for 7*i explicitly in terms of r, and for the remaining variables in terms of 7*i: 



(69) 



v 



l-^(l + K 2 )(lnr + ^ l) + o((^) 2 ) 

C v Ki [l + TZt + (9(7*?)] , (70) 

s + = -iit [i + 2(1 + icDn, + o(n\)] , (71) 

M" 1 = CufKv [1 + §7* a + 0{K\)] , (72) 

as t — > +oo. Here Cm, C v and are constants that depend on the initial conditions. The 

asymptotic expansion for IZ2 can be obtained directly from (|34j) and (|35|). Its leading term 
is 

n 2 = J 1 - f ci n x [\ + cw] , (73) 



which entails the restriction < C v < |. 

We note that the asymptotic expansion for S + also contains a term Cs + e _r , and the 
asymptotic expansion for ip would contain an additive constant C^. Thus all five constants 
of integration for the five-dimensional system (}2"%|) - (j3"4")) have been accounted for, confirming 
that the asymptotic form above describes the general Bianchi VIIo solution with a radiation 
fluid source. We can set one of these constants to unity by using the freedom to shift 
r — > t + C, leaving four essential constants to parametrize the four-dimensional family of 
orbits in the five- dimensional state space. We also note the flawed argument in NHW (see 
page 3126), which should be corrected as above. 

We now give the asymptotic dependence of the physical quantities of primary interest on 
the cosmological clock time t, which is related to r via equation (J2J). In the present situation 
the Raychaudhuri equation reads 

H' = -(l + q)H, (74) 
(see Hewitt et al 2001, equation (A.9)). On account of (|2Tj) . it assumes the form 

#'= -(2 + £;j_ + 7e 2 cos 2V>)# . (75) 
We integrate this equation using (jEHJ), (jZIj) and (fT3~J) . which gives 8 

H= \Ae- 2T [l + 0{l)] . (76) 

8 Note that the integration of the exponentially rapidly oscillating term H2 cos 2ip yields a constant. 
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where A is a constant with dimensions of (time) 1 , and then use (0) to obtain 

t = ie* [1 + 0(1)] . (77) 

Note that on account of (J55j) . ([59]) and (fTTj) . 

n = l + 0(i), (78) 
as t — > +oo. The energy density p can then be determined algebraically through the relation 

n = P /(3H 2 ). 

The desired asymptotic forms, valid for At 3> 1, are as follows: 

H ~ ^- , p ~ -At , u ~ , 7-7 x , (79) 



E 2 



1 



21n(At) 



l + ^/l-fC2 cos 2^ 



10) 



ACUAt) 



(Wi, W 2 ) « - ' A /l - f C v (cos2^,sin2^) , (SI) 



where 



dip Cm A 



dt y/At\n(At) 



i2) 



Observe that the tilt constant C v primarily affects the amplitude of the shear oscillations 
through (J8U|) . while the constant Cm affects how rapidly the shear oscillates through (J82|) . 
and hence the amplitude of the Weyl curvature scalars Wi and W2. 

Although our emphasis in this paper is on the radiation-filled models (7 = |), we also 
give the asymptotic forms as r — > +00 for the models with | < 7 < |, which includes the 
physically important case of dust (7 = 1), in order to make comparisons. The derivation of 
these results is an extension of the analysis given in Appendix B of WHU. 

For any tilted irrotational perfect fluid Bianchi VII cosmology with 1 < 7 < |, 

v = C v e- 2l3T [l + 0(e- bT )] , (83) 
Tlx = C ni e- 2 ^ [1 + 0(e- bT )} , (84) 

TZ 2 = yfa - 7 2 ^ e- 2 ^[l + 0{e^)) , (85) 
S + = [I + 0(e- b n] , (86) 

M" 1 = C M e {1 ~ p)T [l + 0(e- bT )} , (87) 

as t — > +00, where f3 = |(4 — 37), and Cm, C v and C^ are constants that depend on 
the initial conditions, with < C v < and b is a positive constant. We note that 

the asymptotic expansion of S + also contains a term Cs + e~2^ 2_7 - ) ' r , and the asymptotic 
expansion for ip would contain an additive constant C^, giving five constants of integration 
as in the radiation case. 

We now express the asymptotic behaviour of the various physical quantities in terms of 
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clock time t. As in the radiation case we obtain 

H=^e-^\l + 0{e- 2 ^)] , t = ief^ [1 + 0(e- 2 ^ 
37 A 

where A is a constant with dimensions of (time) -1 . 

The desired asymptotic forms, valid for At 3> 1, are as follows: 

H*^, P«g^. (89) 

2(4—3-7) / / \ 

£ 2 » K^-V^ + - 7 2 C 2 cos 2^ J , (90) 

(W 1; W 2 ) « -C 2 M JC1 1 - 7 2 C2(At)^(cos2^sin2^) , (91) 

where 

V>«^(^- (92) 

6 The role of the cosmological constant 

In this section we discuss the influence of a cosmological constant A > on the dynamics 
of tilted irrotational perfect fluid Bianchi VIIo cosmologies. The evolution equations in 
Section El are modified only by the presence of 

«a - ip («» 

in (jm and (JT5J) as follows: 

Q = 1 - T? + - S 2 - S 2 - N 2 _ - iV 2 - tt A (94) 
q = 2(S 2 + + S 2 + S 2 ) + |G; 1 [(3 7 - 2) + (2 - 7 )t; 2 ]fi - Q A , (95) 

and Q A satisfies the evolution equation 

Vl' A = 2(q + l)Vl A . (96) 

The expression for the Weyl curvature components is not affected. 

Wald's theorem (see Wald 1983) makes a strong and general statement about the asymp- 
totic behaviour of spatially homogeneous cosmologies with A > and subject to the stress- 
energy tensor satisfying a physically reasonably energy condition. The result is that the 
models are asymptotic to the de Sitter solution, in the sense that 

lim tt A = 1 , (97) 

T— » + 00 

in our notation. It follows immediately from (J97|) . (|14j) . (J49)) and (|15|) that 

lim (£+, S x , N_, N x ) = 0, (98) 

T— > + 00 
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and hence that 

lim Q = , lim E = , lim q = -1 . (99) 



r— >+oo t— >+oo r^+oo 



Wald's theorem gives no information about the Weyl parameters or about the tilt v. However 
it follows from and JTOJ) that 

lim N + = , (100) 

which implies, on account of (|H3j) - (jfinj) . that 

lim Wi = , lim W 2 = . (101) 



r— >+oo 7" — * -r~ oo 



A significant difference occurs as regards the tilt. It follows from Section III of Lim et al 
2004 that 

'0 if I < 7 < I 



lim v = < 3 4 3 (102) 

^oo if 7 



T— > + 00 



4 

3 ' 



where < < 1, i.e. the value of 7 for a radiation fluid (7 = |) is a bifurcation value for 
the tilt in a model with a cosmological constant. 



7 Discussion 

In this paper we have derived the late-time dynamical behaviour of tilted Bianchi VIIo 
cosmologies with an irrotational radiation fluid as source, and have given the asymptotic 
form of the general solution as t —>■ +00. 

The main result is Theorem^ which establishes rigorously the asymptotic behaviour of 
the cosmological models. 9 Theorem which is an immediate consequence of Theorem Q 
describes the degree to which the models isotropize. These theorems generalize the results 
of NHW for the corresponding non-tilted models. The derivation of the detailed asymptotic 
form of the solutions as t — > +00, given in Section is not completely rigorous, since it 
relies on the analysis of a so-called truncated system of differential equations, as in NHW. 

This class of models displays several noteworthy features. Firstly, the models appear 
to approximate the flat FL model at late times, since the Hubble-normalized shear and 
the tilt tend to zero, and the density parameter tends to one. The Hubble-normalized 
Weyl curvature scalars Wi and W2 diverge, however, indicating that physically significant 
anisotropy remains. In effect, the models isotropize as regards shear and tilt, but not as 
regards Weyl curvature. The scalars Wi, W2, which are defined by equations (|5Uj) and (pTT|). 
diverge and oscillate increasingly rapidly in an out-of-phase manner (see equation 1)81))). The 
scalar W, defined by 

W 2 = W 2 + W 2 2 , (103) 

then satisfies 

lim W = +00 , (104) 

as follows from (|U3*|) and ([64)1 . Thus, in terms of Hubble-normalized scalars, the future 
9 The assumption IZ2 > of Theorem ^ is supported numerically (see paragraph before (|37|l V 
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asymptotic state (t — > +00) for a radiation-filled universe is given by 

E ^ , v^O , Q -»• 1 , W ^ 00 . (105) 

Secondly, the isotropization as regards shear and tilt occurs at a slow rate, i.e. E 2 and v 
tend to zero logarithmically in terms of clock time t (proportional to , see equations 
dZHJ) and (JBOJ))- 

The dust-filled (7 = 1) universes differ significantly. The Weyl curvature scalar W 
remains bounded and approaches a finite value Woo that depends on the initial conditions 
(see equation (|§Tj) with 7 = 1). Thus, in terms of Hubble-normalized scalars, the future 
asymptotic state for a dust-filled universe is given by 

E -»• , v^O , -»• 1 , W^Woo. (106) 

The isotropization as regards shear and tilt occurs at a faster rate, i.e. E 2 and v tend 
to zero like (At)~ 2 ^ 3 (see equation (|90|) with 7 = 1). We note in passing that a universe 
described by (|106|) . i.e. one that is highly isotropic as regards the shear but not as regards 
the Weyl curvature, is compatible with the observed highly isotropic temperature of the 
cosmic microwave background radiation (see WHU, Section 4, for further discussion of this 
point). 10 

As shown in Section 03 the dominance of the Weyl curvature at late times is curbed by 
the presence of a positive cosmological constant (see equation (JSHJl), i-e. the Weyl scalar 
W tends to zero. Thus, in terms of Hubble-normalized scalars, the future asymptotic state 
(t — > +00) for a universe with a cosmological constant and radiation or dust, is 

E ^ , v^ Voo , fi A 1 , ^^0, W^O, (107) 

with v 00 = for dust and < v 00 < 1 for radiation. The asymptotic states ()105|) or ()106J) 
with A = are still of relevance for models with A > since they can act as intermediate 
states, i.e. a model with A > can be approximated arbitrarily closely by the states A 10 5}) or 
MUb]) over a finite time interval. 

In conclusion, we note that the most striking feature of the late-time dynamics of Bianchi 
VIIo cosmologies with a radiation fluid is the so-called Weyl curvature dominance: if the 
cosmological constant is zero (A = 0) then the Hubble-normalized Weyl scalar W diverges 
as t — > +00, while if A > 0, the scalar W can become arbitrarily large during a finite 
time interval. The phenomenon was first encountered for non-tilted Bianchi VII models 
(see WHU), and in this paper we have shown that it persists in the presence of tilt, for an 
irrotational fluid. The same phenomenon has recently been found in Bianchi VIIo models 
with a tilted and rotating fluid (Hervik et al 2006). Weyl curvature dominance also occurs 
in non-tilted Bianchi VIII cosmologies (Horwood et al 2003), and has recently been found 
in tilted Bianchi VIII cosmologies (Hervik & Lim 2006). 

10 Dominance of the Weyl curvature may lead to other observable effects, such as image distortion (see, 
e.g., Chrobok & Perlick 2001. 
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A Proof of Theorem [T] 

In the proof we need the following standard result and a second lemma. 

Lemma 1 Let /(r) be a non-negative real-valued and Lipschitz continuous function on the 
interval r < r < +00. If 

POO 

f(r) dr is finite 



TO 



then lim r ^ +00 /(r) = 0. 

Lemma 2 For solutions of equations i28\) - fifty) , the following limits cannot both hold: 

lim H = , lim E, = a , < a < 1. (108) 



T^+oo T^+oo 



Proof. Suppose that both hold. Then ()35|) implies that lim T ^ +00 7?.i = 1 — a 2 . It follows 
from ()38|) and (pTT|) that lim T _» +00 S + = a, lim T _^ +00 f2 = 0. Then (|3*7j) and (}3*Hj) imply that 
lim T ^ +00 T,' + = —(1 + a)(l — a 2 ) 7^ 0, which contradicts lim T ^ +00 S + = a. □ 

We begin by proving that 

lim v = , (109) 

T— > + 00 

using the evolution equations for S + and t> in the form (|2*8|) and which we repeat for 
convenience: 



3-v 2 



3 + v 2 ' 

E + (l-v 2 )v. (Ill) 



E' + = -(2 - g)S + - (iV! + iV 2 ) - ^-^ 2 ( I 1 0) 

6 



Equation (|110|) implies that if S + (r ) = for some r , then S' + (r ) < 0, and S + < for 
r > T . We conclude that 

the sign of E + is fixed eventually. (112) 
Hence v is monotone eventually and since it is also bounded, 

lim v = L , (113) 

where < L < 1. Hence 

00 

v dr is finite. (114) 

TO 
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It follows from ()110|) - (jlllj) that v" is bounded, which implies that v' is Lipschitz continuous. 
We thus conclude from Lemma Q that 



lim v' = . (115) 
We want to show that L = 0. To do this we assume the contrary, namely that < L < 1. 
Case 1: L = 1 

Since v < 1 and is monotone eventually then t> is eventually increasing. Hence £ + > 
eventually by (|TTTj) and ()112j) . It follows from ()110|) that £' + < 0, hence £+ is monotone 
decreasing eventually, and so 

lim £+ = a , < a < 1 . (116) 

It follows that 



T^OO 



TO 



E' + dr is finite. (117) 
Since £+ > and L ^ it follows from (jll(J|) that 

< i^f) ^ C, (" S V) f ^ r > r , (118) 

where C is a positive constant, and tq is sufficiently large. We conclude that Q dr is 
finite. Since Q' is bounded by f2 is Lipschitz continuous, and it follows from Lemma [T] 
that 

lim Q = . (119) 

By Lemma El we have a contradiction. 
Case 2: < L < 1 

It follows from ffTTTft . (fTT3|l and fTTKjl that lim T ^ 00 S + = 0, and by integrating (fTTTJl it 
follows that £ + dr is finite. Since L > 0, it follows from (jllOj) that 

<-[£; + (2 - 9 )£ + ]f^±^ <-C[£' + + (2- ? )£ + ] for r>r , (120) 

where C is a positive constant, and To is sufficiently large. Since 2 — q is bounded, we conclude 
that J" f2 dr is finite. The rest of the proof is the same as case 1. 

To summarize, we have established that ()109|) holds. The next step is to show that 
lim T ^ +QO Q exists. Extending the method used in Appendix A of NHW we integrate to 
obtain 

i, n . r4£,v 2 



l n?r = / £^rfr+ / T -^dr+ / MK^ dr , (121) 

z » '0 Jt </t d + U -/t 

where = ^(tq) and the initial time is chosen such that M < 1 for all r > r . 



16 



Following Section 3 of NHW, we introduce an auxiliary variable Z according to 



Vn 



(122) 



whose evolution equation is 
Z' 



-1 -2(1 + E+) — - cos 2^ - 3ft 2 M sin 2?/> ^S, 

Tci 3 + v 2 ' 



(123) 



as follows from (|2*5j) . (|3*T|) and (|3l)J). We modify Z according to 

Z 



l-|M(l + S + )|a s in2^ 



and obtain 



Z' 



3 + v 2 



Since v and M tend to zero as r — > +oo, equation ()125|) implies that 

z = e>(e ( - 1+5)r ) 

where 5 can be chosen arbitrarily small. Equation ()122|) and ()124j) then imply that 

POO 

MTl x = 0(e { ~ 1+s)T ) and hence / MTLxB^ dr is finite. 



From (jllljl we obtain 



4E,t> 2 



2(3 - f 2 )ff ' 



3 + w 2 3(3 + w 2 )(l-t; 2 ) 



(124) 



(125) 



(126) 



(127) 



128) 



where f(v) = | In j^tt, and hence the second integral in ()121j) is given by 



(129) 



where v = v(t ). It follows from (jlU9j) that 

4E+^ 2 
3 + w 2 



dr is finite. 



(130) 
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We have thus shown that the second and third integral terms in ()121j) are bounded at late 
times. Since Q is bounded above, it follows from (|121jl that 



Kr) :-- 



E+ dr is bounded above for all r > r . 



(131) 



Since the integrand El is non-negative, the function I(t) is monotonically increasing, which 
implies that lim T ^ +00 I(t) exists and is finite. It follows from (j!21j) . (J127)) and ()13())1 that 
limT-^+oo Q exists, and hence from (|4*T|) that lim r _^ +00 Q exists, i.e. there is a constant L, 
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satisfying < L < 1 such that 

lim tt = L . (132) 



T— > + 00 



The next step is to show that E + and TZi tend to zero and that L = 1. Since S + and TZi 
are bounded, it follows from ()42|) that (£+)' = 2E + E' + is uniformly bounded for r > ro. This 
implies that T? + is Lipschitz continuous for r > To- It thus follows from (|131|) and Lemma ^ 
that 

lim S + = , (133) 



T— > + 0O 



and hence, from (J58|) and (J37J), that 

lim £+ = . (134) 

It now follows from ([35 )1 . (|59 j) . ffrTT fl . (JHU) and (O that 

lim ^! = lim 7^ = 1 - L . (135) 

The limits (|3ZD, fTTEty . ifTEty . fT33j) and equation (H2J imply that 

lim E' = L - 1 , (136) 

which contradicts ()133j) unless L = 1. Using ()135|) we obtain 

lim ^i = 0. (137) 

T— > + 0O 

It follows from flHU) and (ITUHj) that 

lim TZ 2 = . (138) 

T— > + 0O 

Next, we will show that M 2 /Ki tends to zero. We note that (pE3j) - p4]l give 

fM 2 \' - M 2 

' —J = "2 [1 + 3S + + 2£* + M(Bm + B-nJ] ^- . (139) 



It then follows from (|133|) and (|37j) that 



M 2 

lim = , (140) 



and hence that 



M 2 1Z\ 
lim — — = , or equivalently lim — — = +oo . (141) 

T^ + OO TCI ~ T-+ + 00 

It remains to show that 

lim^ = +oo. (142) 
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Equation ()33|) can be written in the form 

V = 2E+V , (143) 

where 

V = — ^ t • (144) 

(1 — V Z ) 3 

Let 

Y=Z-. (145) 



It follows from (|38j) . (j43 jl and (JT43J) that 

V = -2(E^ + MB)Y, (146) 
where 5 is bounded as r — > +oo. This equation can be integrated to yield 

Y = Y e Jt o + e Jt o . (147) 

We know that 

lim / Si dr is finite. (148) 

T— > + 00 



TO 



In addition equations (jHJ), and (|133|) imply that 

M = C(e ( - 1+(5)r ), as r -> +oo, (149) 

where 5 is an arbitrarily small positive number. It thus follows from (I145|) and (|147|) - (jl49j) 

that 

Urn ¥- = C, (150) 
which, with equations (|143|) . and (J3~7j) . implies that 



lira ^ = (7 . (151) 



We now write (J3*%j) in the form 



It follows from ()134|) . ()137f) and (J33j) that lim T ^ +00 = 1. Equation (j!51|) thus implies 

^(l) 2 = 1 - fc2 ' (153) 

which is non-zero in general. Finally, since 

*. /'R.A/Ki. (154) 



m 2 V^iy v^ 2 
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it follows from (|141|) that 

m = +°°- ( 155 ) 

□ 



B Derivation of decay rates 

In this Appendix we use centre manifold theory (see for example, Carr 1981, and NHW, Ap- 
pendix B, for a brief summary) to derive the asymptotic decay rates of solutions of equations 
(foT)j) - (p)%j) as r — > +00. The equilibrium point (E + ,R,v) = (0,0,0) is asymptotically stable, 
but non-hyperbolic, having two zero eigenvalues (those associated with R and v). There is 
thus a two-dimensional centre manifold of the form 



h(R,v) 



(156) 



Using the method summarized in NHW, Appendix B, we find the following first order ap- 
proximation to the centre manifold: 



S + = — R 2 — \y 2 + {higher order terms} 



(157) 



To determine the asymptotic behaviour on the centre manifold, we observe that (|68j) can be 
written as 

v 



V = 2£ , V 



V 



which gives 



(1 _ #2)1/3 

2±Us) ds 



V. J TO 



V(r) = V(r )exp 
Since lim r _» +00 V{r) = 0, we obtain 

lim / S+(s) ds = —00 . 

Now observe that (|67|) gives 

R(t) = R{t ) exp 



(158) 
(159) 

(160) 









/ S + ds 


exp 


/ t\ ds 


-J TO 




J TQ 



(161) 



where the second exponential term is greater than unity, since the integrand is positive. 
Hence 



V _ V{t ) 

r 2 R(r y 



exp 



SI ds 



TO 



and 



lim — = C 



(162) 
(163) 
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i.e. the limit exists, but could be zero. In other words, v = 0(R 2 ). So the first order 
approximation of the centre manifold (|157|) becomes 



E + = — R 2 + {higher order terms} . 
We can use this information to integrate (|67j) . obtaining 

R 2 = h {higher order terms} , 

2r 



and so 



which leads to 



At 2 



±lds = C- — + o(- 
Ar \ r 



2 / ' 



C - \R 2 + o{R 2 



where C is a positive constant. Hence from (j!58|) and (j!62|) we obtain 

v = C V {R 2 + R 4 ) + o(i? 4 ) , 



(164) 

(165) 

(166) 
(167) 

(168) 



where C v is an arbitrary constant. We can now reapply the centre manifold procedure to 
amend ()157j) by including higher order terms that are comparable to 



v 2 : 



E + = —R 2 — 2i? 4 — |t> 2 + {higher order terms} . 
We then solve (ffiTj) to obtain 



R 



1 



l + IC 2 )(\nr + C ni : 
At 



or 



(169) 



(170) 



and substitute into ()168j) and ()157|) to give the asymptotic form of v and E+. The error 
bound o{t~ 1 ) can be improved to O ( ( — ) 2 ) , by increasing the accuracy of the centre manifold 
approximation (see NHW, page 3132). 
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